Genomic insights into the Montseny brook newt (Calotriton arnoldi), a Critically Endangered glacial relict.

The Montseny brook newt (Calotriton arnoldi), considered the most endangered amphibian in Europe, is a relict salamandrid species endemic to a small massif located in northeastern Spain. Although conservation efforts should always be guided by genomic studies, those are yet scarce among urodeles, hampered by the extreme sizes of their genomes. Here, we present the third available genome assembly for the order Caudata, and the first genomic study of the species and its sister taxon, the Pyrenean brook newt (Calotriton asper), combining whole-genome and ddRADseq data. Our results reveal significant demographic oscillations which accurately mirrored Europe's climatic history. Although severe bottlenecks have led to depauperate genomic diversity and long runs of homozygosity along a gigantic genome, inbreeding might have been avoided by assortative mating strategies. Other life history traits, however, seem to have been less advantageous, and the lack of land dispersal has driven to exceptional levels of population fragmentation.

Hidden in the sand: Phylogenomics unravel an unexpected evolutionary history for the desert-adapted vipers of the genus Cerastes.

The desert vipers of the genus Cerastes are a small clade of medically important venomous snakes within the family Viperidae. According to published morphological and molecular studies, the group is comprised by four species: two morphologically similar and phylogenetically sister taxa, the African horned viper (Cerastes cerastes) and the Arabian horned viper (Cerastes gasperettii); a more distantly related species, the Saharan sand viper (Cerastes vipera), and the enigmatic Böhme's sand viper (Cerastes boehmei), only known from a single specimen in captivity allegedly captured in Central Tunisia. In this study, we sequenced one mitochondrial marker (COI) as well as genome-wide data (ddRAD sequencing) from 28 and 41 samples, respectively, covering the entire distribution range of the genus to explore the population genomics, phylogenomic relationships and introgression patterns within the genus Cerastes. Additionally, and to provide insights into the mode of diversification of the group, we carried out niche overlap analyses considering climatic and habitat variables. Both nuclear phylogenomic reconstructions and population structure analyses have unveiled an unexpected evolutionary history for the genus Cerastes, which sharply contradicts the morphological similarities and previously published mitochondrial approaches. Cerastes cerastes and C. vipera are recovered as sister taxa whilst C. gasperettii is a sister taxon to the clade formed by these two species. We found a relatively high niche overlap (OI > 0.7) in both climatic and habitat variables between C. cerastes and C. vipera, contradicting a potential scenario of sympatric speciation. These results are in line with the introgression found between the northwestern African populations of C. cerastes and C. vipera. Finally, our genomic data confirms the existence of a lineage of C. cerastes in Arabia. All these results highlight the importance of genome-wide data over few genetic markers to study the evolutionary history of species.

Evolution along allometric lines of least resistance: morphological differentiation in Pristurus geckos.

Species living in distinct habitats often experience unique ecological selective pressures, which can drive phenotypic divergence. However, how ecophenotypic patterns are affected by allometric trends and trait integration levels is less well understood. Here we evaluate the role of allometry in shaping body size and body form diversity in Pristurus geckos utilizing differing habitats. We found that patterns of allometry and integration in body form were distinct in species with different habitat preferences, with ground-dwelling Pristurus displaying the most divergent allometric trend and high levels of integration. There was also strong concordance between intraspecific allometry across individuals and evolutionary allometry among species, revealing that differences in body form among individuals were predictive of evolutionary changes across the phylogeny at macroevolutionary scales. This suggested that phenotypic evolution occurred along allometric lines of least resistance, with allometric trajectories imposing a strong influence on the magnitude and direction of size and shape changes across the phylogeny. When viewed in phylomorphospace, the largest rock-dwelling species were most similar to the smallest ground-dwelling species, and vice versa. Thus, in Pristurus, phenotypic evolution along the differing habitat-based allometric trajectories resulted in similar body forms at differing body sizes in distinct ecological habitats.

Genomics reveals introgression and purging of deleterious mutations in the Arabian leopard (Panthera pardus nimr).

In endangered species, low-genetic variation and inbreeding result from recent population declines. Genetic screenings in endangered populations help to assess their vulnerability to extinction and to create informed management actions toward their conservation efforts. The leopard, Panthera pardus, is a highly generalist predator with currently eight different subspecies. Yet, genomic data are still lacking for the Critically Endangered Arabian leopard (P. p. nimr). Here, we sequenced the whole genome of two Arabian leopards and assembled the most complete genomic dataset for leopards to date. Our phylogenomic analyses show that leopards are divided into two deeply divergent clades: the African and the Asian. Conservation genomic analyses indicate a prolonged population decline, which has led to an increase in inbreeding and runs of homozygosity, with consequent purging of deleterious mutations in both Arabian individuals. Our study represents the first attempt to genetically inform captive breeding programmes for this Critically Endangered subspecies.

Species boundaries to the limit: Integrating species delimitation methods is critical to avoid taxonomic inflation in the case of the Hajar banded ground gecko (Trachydactylus hajarensis).

With the advent of molecular phylogenetics and the implementation of Multispecies Coalescent-based (MSC) species delimitation methods (SDM), the number of taxonomic studies unveiling and describing cryptic sibling species has greatly increased. However, speciation between early divergent lineages is often defined without evaluating population structure or gene flow, which can lead to false claims of species status and, subsequently, taxonomic inflation. In this study, we focus on the intriguing case of the Arabian gecko Trachydactylus hajarensis (Squamata: Gekkonidae). We generated mitochondrial data (12S rDNA) and genome-wide SNP data (ddRADseq) for 52 specimens to determine phylogenomic relationships, population structure and genetic diversity within this species. Then, we applied a set of different SDMs to evaluate several competing species hypotheses through the MSC. Results show that T. hajarensis is comprised by three well-defined population lineages, two of them in the Hajar Mountains of eastern Arabia, and one in Masirah Island, on the southeastern coast of Oman. Strong mito-nuclear discordances were found between populations inhabiting the Hajar Mountains, but we did not find evidence of current gene flow between them. Surprisingly, discordances in species tree topology were found when different downsampled datasets were used, and especially when linking population sizes, a commonly implemented feature in species tree reconstruction with genomic data. Different SDMs yielded different results, supporting from four species within the group, to T. hajarensis being a single species. With such contrasting results we suggest caution before splitting T. hajarensis. Overall, this study highlights the importance of sample and prior choice and the integration of several SDMs to not incur into taxonomic inflation, providing a set of already available tools to assess population structure, genetic diversity, and SDMs before describing new species.

Reconstructing Squamate Biogeography in Afro-Arabia Reveals the Influence of a Complex and Dynamic Geologic Past.

The geographic distribution of biodiversity is central to understanding evolutionary biology. Paleogeographic and paleoclimatic histories often help to explain how biogeographic patterns unfold through time. However, such patterns are also influenced by a variety of other factors, such as lineage diversification, that may affect the probability of certain types of biogeographic events. The complex and well-known geologic and climatic history of Afro-Arabia, together with the extensive research on reptile systematics in the region, makes Afro-Arabian squamate communities an ideal system to investigate biogeographic patterns and their drivers. Here, we reconstruct the phylogenetic relationships and the ancestral geographic distributions of several Afro-Arabian reptile clades (totaling 430 species) to estimate the number of dispersal, vicariance and range contraction events. We then compare the observed biogeographic history to a distribution of simulated biogeographic events based on the empirical phylogeny and the best-fit model. This allows us to identify periods in the past where the observed biogeographic history was likely shaped by forces beyond the ones included in the model. We find an increase in vicariance following the Oligocene, most likely caused by the fragmentation of the Afro-Arabian plate. In contrast, we did not find differences between observed and expected dispersal and range contraction levels. This is consistent with diversification enhanced by environmental processes and with the establishment of a dispersal corridor connecting Africa, Arabia and Eurasia since the middle Miocene. Finally, here we show that our novel approach is useful to pinpoint events in the evolutionary history of lineages that might reflect external forces not predicted by the underlying biogeographic model. [Dispersal; diversification; model adequacy; paleogeography; reptiles; simulations; vicariance.].

Ecological specialization, rather than the island effect, explains morphological diversification in an ancient radiation of geckos.

Island colonists are often assumed to experience higher levels of phenotypic diversification than continental taxa. However, empirical evidence has uncovered exceptions to this 'island effect'. Here, we tested this pattern using the geckos of the genus Pristurus from continental Arabia and Africa and the Socotra Archipelago. Using a recently published phylogeny and an extensive morphological dataset, we explore the differences in phenotypic evolution between Socotran and continental taxa. Moreover, we reconstructed ancestral habitat occupancy to examine if ecological specialization is correlated with morphological change, comparing phenotypic disparity and trait evolution between habitats. We found a heterogeneous outcome of island colonization. Namely, only one of the three colonization events resulted in a body size increase. However, in general, Socotran species do not present higher levels or rates of morphological diversification than continental groups. Instead, habitat specialization explains better the body size and shape evolution in Pristurus. Particularly, the colonization of ground habitats appears as the main driver of morphological change, producing the highest disparity and evolutionary rates. Additionally, arboreal species show very similar body size and head proportions. These results reveal a determinant role of ecological mechanisms in morphological evolution and corroborate the complexity of ecomorphological dynamics in continent-island systems.

Diversity, multifaceted evolution, and facultative saprotrophism in the European Batrachochytrium salamandrivorans epidemic.

While emerging fungi threaten global biodiversity, the paucity of fungal genome assemblies impedes thoroughly characterizing epidemics and developing effective mitigation strategies. Here, we generate de novo genomic assemblies for six outbreaks of the emerging pathogen Batrachochytrium salamandrivorans (Bsal). We reveal the European epidemic currently damaging amphibian populations to comprise multiple, highly divergent lineages demonstrating isolate-specific adaptations and metabolic capacities. In particular, we show extensive gene family expansions and acquisitions, through a variety of evolutionary mechanisms, and an isolate-specific saprotrophic lifecycle. This finding both explains the chytrid's ability to divorce transmission from host density, producing Bsal's enigmatic host population declines, and is a key consideration in developing successful mitigation measures.

Combined Molecular and Elemental Mass Spectrometry Approaches for Absolute Quantification of Proteomes: Application to the Venomics Characterization of the Two Species of Desert Black Cobras, Walterinnesia aegyptia and Walterinnesia morgani.

We report a novel hybrid, molecular and elemental mass spectrometry (MS) setup for the absolute quantification of snake venom proteomes shown here for two desert black cobra species within the genus Walterinnesia, Walterinnesia aegyptia and Walterinnesia morgani. The experimental design includes the decomplexation of the venom samples by reverse-phase chromatography independently coupled to four mass spectrometry systems: the combined bottom-up and top-down molecular MS for protein identification and a parallel reverse-phase microbore high-performance liquid chromatograph (RP-µHPLC) on-line to inductively coupled plasma (ICP-MS/MS) elemental mass spectrometry and electrospray ionization quadrupole time-of-flight mass spectrometry (ESI-QToF MS). This allows to continuously record the absolute sulfur concentration throughout the chromatogram and assign it to the parent venom proteins separated in the RP-µHPLC-ESI-QToF parallel run via mass profiling. The results provide a locus-resolved and quantitative insight into the three desert black cobra venom proteome samples. They also validate the units of measure of our snake venomics strategy for the relative quantification of snake venom proteomes as % of total venom peptide bonds as a proxy for the % by weight of the venom toxins/toxin families. In a more general context, our work may pave the way for broader applications of hybrid elemental/molecular MS setups in diverse areas of proteomics.

Small and overlooked: Phylogeny of the genus Trigonodactylus (Squamata: Gekkonidae), with the first record of Trigonodactylus arabicus from Jordan.

Geckos of the genus Trigonodactylus are widely distributed in the sand deserts of the Arabian Peninsula. Three species of this genus are currently recognized, with a fourth one, Stenodactylus pulcher, which placement within Trigonodactylus has been tentatively suggested, but not yet confirmed. We present a phylogenetic analysis of the genus Trigonodactylus with new specimens collected in central Saudi Arabia and southern Jordan. New genetic data has been generated from three mitochondrial markers to investigate the phylogenetic relationships of all species of the genus and to assess the putative generic assignment of S. pulcher. Our results confirm that S. pulcher indeed belongs within Trigonodactylus, branching as a sister lineage to all other species of the genus. The new samples cluster within Trigonodactylus arabicus, thus confirming the genetic homogeneity of the species across its large and seemingly inhospitable range. The new specimen collected in southern Jordan represents the first record for the country and a considerable range extension to the northwest from all previously reported localities. Our findings and discovery of a new species for Jordan highlight the need of more field surveys to be carried out in the underexplored parts of Jordan and northern Saudi Arabia, as these places still hold a potential for new discoveries and are crucial for understating the biogeography of the Arabian herpetofauna.

Contributions to ophiology in Zootaxa 20012020: patterns and trends.

There are currently 3,900 recognized, extant snake species belonging to 529 genera globally (Uetz et al. 2021; this study), making snakes one of the most diverse major groups of squamates. Of the 665 currently recognized species that were described between 2001 and 2020 (a ~17% increase in total species), ~34% of these (226 species) were described in Zootaxa. This number does not include species resurrected from synonymy. The other ~66% (439) species were described in 105 other journals, bulletins or books (Fig.1a). Overall, the number of new snake species described every year is gradually increasing, and 40% of the new species described since 2011 were published in Zootaxa. Following Zootaxa, the second ranked journal, with 37 described species since 2001, is Herpetologica (Fig. 2). Anecdotally, the choice of Zootaxa as a publication outlet for new species descriptions by most authors is based on speed of publication post-acceptance, publication free of charge, relatively unconstrained number of papers published per year, relatively unconstrained manuscript length, expert section editors and reviewers, and consolidated scientometric parameters.

The demise of a wonder: Evolutionary history and conservation assessments of the Wonder Gecko Teratoscincus keyserlingii (Gekkota, Sphaerodactylidae) in Arabia.

Effective biodiversity conservation planning starts with genetic characterization within and among focal populations, in order to understand the likely impact of threats for ensuring the long-term viability of a species. The Wonder Gecko, Teratoscincus keyserlingii, is one of nine members of the genus. This species is distributed in Iran, Afghanistan, and Pakistan, with a small isolated population in the United Arab Emirates (UAE), where it is classified nationally as Critically Endangered. Within its Arabian range, anthropogenic activity is directly linked to the species' decline, with highly localised and severely fragmented populations. Here we describe the evolutionary history of Teratoscincus, by reconstructing its phylogenetic relationships and estimating its divergence times and ancestral biogeography. For conservation implications of T. keyserlingii we evaluate the genetic structure of the Arabian population using genomic data. This study supports the monophyly of most species and reveals considerable intraspecific variability in T. microlepis and T. keyserlingii, which necessitate broad systematic revisions. The UAE population of T. keyserlingii likely arrived from southern Iran during the Pleistocene and no internal structure was recovered within, implying a single population status. Regional conservation of T. keyserlingii requires improved land management and natural habitat restoration in the species' present distribution, and expansion of current protected areas, or establishment of new areas with suitable habitat for the species, mostly in northern Abu Dhabi Emirate.

Systematics, biogeography and evolution of the Saharo-Arabian naked-toed geckos genus Tropiocolotes.

Plate tectonics constitute one of the main mechanisms of biological diversification on Earth, often being associated with cladogenetic events at different phylogenetic levels, as well as with exchange of faunas and floras across previously isolated biogeographic regions. North Africa and Arabia share a complex geological history that dates back to the break-up of the Arabian plate from the African plate ~30-25 Mya, followed by various geological events, such as the formation of the Red Sea or the connection between the African, Arabian and Eurasian plates. Species with Saharo-Arabian distributions have shown a close association between their evolutionary history and these geological events. In this study, we investigate the systematics, biogeography and evolution of the genus Tropiocolotes, a group of small ground-dwelling geckos, comprised by 12 species distributed from the Atlantic coast of North Africa to southwestern Iran. Species delimitation analyses uncovered the existence of high levels of undescribed diversity, with forms here considered at the species level including Tropiocolotes tripolitanus (Mauritania and southern Morocco), T. nattereri (southern Israel) and T. scorteccii (Yemen and Oman). Phylogenetic and biogeographic analyses recovered two main clades, an exclusively African clade and a Saharo-Arabian clade, that split ~25 Mya following the vicariant event mediated by the separation of the Arabian and African plates. The complex geological activity around the Red Sea is associated with the diversification within the Saharo-Arabian clade, including the colonization of North Africa from a second Tropiocolotes group. Results also provide new insights into the geographic distribution of Tropiocolotes nubicus, previously considered as exclusively associated to the Nile River valley, extending its known distribution further west, up to the Central Mountains of the Sahara. Accordingly, the Nile River seems to act as a major biogeographic barrier, separating Tropiocolotes nubicus and T. steudneri in their western and eastern margins, respectively.

Parasites in a hotspot: diversity and specificity patterns of apicomplexans infecting reptiles from the Socotra Archipelago.

Although parasites represent a major component of biodiversity, they remain poorly assessed, especially in remote regions. In this study, we screened 461 reptiles from Socotra, the largest and most biologically diverse archipelago in Arabia. Using 18S rRNA primers, we detected various apicomplexan parasites, namely haemogregarines, sarcocystids and eimeriids. Haemogregarines were the most common and genetically diverse, followed by sarcocystids (genus Sarcocystis) and eimeriids (genera Isospora and Lankesterella). All were related to parasites of other reptiles, including species from Arabia, Northern Africa and Asia. Like their 29 endemic reptile hosts, almost all Socotran parasites presented high genetic divergence and ecological differences from those found elsewhere, and probably represent undescribed endemic species. Among hosts, skinks were the most parasitized, which contrasted with similar studies from other areas, probably due to their more generalist diet and habitat use. As expected due to its high species richness, geckos harboured the highest parasite diversity in the archipelago. Parasite diversity also seemed to be correlated to island size, as the largest island harboured most haplotypes. This study emphasizes the importance of screening parasites in wild hosts from remote regions and of considering host ecology to understand disease transmission across taxa.

Swimming through the sands of the Sahara and Arabian deserts: Phylogeny of sandfish skinks (Scincidae, Scincus) reveals a recent and rapid diversification.

Large parts of the Sahara Desert and Arabia are covered by sand seas and sand dunes, which are inhabited by specialized animal communities. For example, many lizards have developed adaptations to life in loose sand, including sand-swimming behavior. The best-known sand swimmers of the Saharo-Arabia are the sandfish skinks (genus Scincus). Although there are currently only four Scincus species recognized, their phylogenetic relationships have not yet been addressed in detail. We use eight genetic markers (three mitochondrial, five nuclear) and a complete sampling of species to infer the relationships within the genus. We employ multiple phylogenetic approaches to reconstruct the evolutionary history of these skinks and to assess the level of reticulation at the onset of their radiation. Our results indicate the presence of five strongly supported species-level lineages, four represented by the currently recognized species and the fifth by S. scincus conirostris, which does not form a clade with S. scincus. Based on these results we elevate the Iranian and northern Arabian S. conirostris to the species level. The two Saharan species, S. albifasciatus and S. scincus, are sister in all analyses. Deeper relationships within the genus, however, remained largely unresolved despite the extensive genetic data set. This basal polytomy, together with the fact that we detected no sign of hybridization in the history of the genus, indicates that the diversification of the five Scincus species was rapid, burst-like, and not followed by secondary hybridization events. Divergence time estimations show a Middle Pliocene crown radiation of the genus (3.3 Mya). We hypothesize that the aridification of the Saharo-Arabia that began in the Late Miocene triggered the initial diversification of Scincus, and that the subsequent expansion of sand deserts enabled their dispersal over the large Saharan and Arabian range. We discuss the evolution of body form in sand swimming lizards and ponder how Scincus retained their fully limbed morphology despite being sand swimmers that are typically limbless.

Presence of low virulence chytrid fungi could protect European amphibians from more deadly strains.

Wildlife diseases are contributing to the current Earth's sixth mass extinction; one disease, chytridiomycosis, has caused mass amphibian die-offs. While global spread of a hypervirulent lineage of the fungus Batrachochytrium dendrobatidis (BdGPL) causes unprecedented loss of vertebrate diversity by decimating amphibian populations, its impact on amphibian communities is highly variable across regions. Here, we combine field data with in vitro and in vivo trials that demonstrate the presence of a markedly diverse variety of low virulence isolates of BdGPL in northern European amphibian communities. Pre-exposure to some of these low virulence isolates protects against disease following subsequent exposure to highly virulent BdGPL in midwife toads (Alytes obstetricans) and alters infection dynamics of its sister species B. salamandrivorans in newts (Triturus marmoratus), but not in salamanders (Salamandra salamandra). The key role of pathogen virulence in the complex host-pathogen-environment interaction supports efforts to limit pathogen pollution in a globalized world.

Cryptic diversity and non-adaptive radiation of montane New Guinea skinks (Papuascincus; Scincidae).

New Guinea, the world's largest and highest tropical island, has a rich but poorly known biota. Papuascincus is a genus of skinks endemic to New Guinea's mountain regions, comprising two wide-ranging species and two species known only from their type series. The phylogeny of the genus has never been examined and the relationships among its species - as well as between it and closely related taxa - are hitherto unknown. We performed the first large-scale molecular-phylogenetic study of Papuascincus, including sampling across the genus' range in Papua New Guinea. We sequenced three mitochondrial and two nuclear markers from 65 specimens of Papuascincus and reconstructed their phylogenetic relationships. We also performed species-delimitation analyses, estimated divergence times and ancestral biogeography, and examined body-size evolution within the genus. Papuascincus was strongly supported as monophyletic. It began radiating during the mid-Miocene in the area now comprising the Central Cordillera of New Guinea, then dispersed eastward colonising the Papuan Peninsula. We found evidence of extensive cryptic diversity within the genus, with between nine and 20 supported genetic lineages. These were estimated using three methods of species delimitation and predominantly occur in allopatry. Distribution and body-size divergence patterns indicated that character displacement in size took place during the evolutionary history of Papuascincus. We conclude that the genus requires comprehensive taxonomic revision and likely represents a species-rich lineage of montane skinks.

Using multilocus approach to uncover cryptic diversity within Pseudotrapelus lizards from Saudi Arabia.

The agamid Pseudotrapelus lizards inhabit the mountainous areas of the Arabian Peninsula and eastern North Africa. Currently six Pseudotrapelus species are recognised, though diagnostic morphological characters are still lacking, creating great difficulty in describing new species. Recently, two specimens of Pseudotrapelus were collected from the vicinity of Riyadh in central Saudi Arabia, an area that was not sampled in previous phylogenetic studies. In here we used both mitochondrial and nuclear data to investigate the phylogenetic position of the new samples, and assess their phylogenetic relationships with the other recognised species of Pseudotrapelus from across the distribution range of the genus. We used a multilocus approach of haplotype networks, concatenated datasets and species trees, performed mitochondrial and nuclear species delimitation analyses, and estimated divergence times. In general, our results support previous molecular studies and uncover the presence of cryptic diversity within Pseudotrapelus. The phylogenetic structure of the genus is of two major clades and within them seven distinct, delimited phylogenetic groups belonging to the six recognised species and the seventh to the individuals from Riyadh. The Riyadh specimens were distinct in all analyses performed. We suggest that the new specimens from the Riyadh area are a distinct lineage, forming a clade with their phylogenetic relatives, P. sinaitus and P. chlodnickii. The clade formed by these three species diverged during the Late Miocene around 6.4 Ma, with cladogenesis possibly facilitated by vicariance and isolation caused due to climatic fluctuations and the progression of sandy areas. Our results suggest further morphological research is necessary to revise the taxonomic status of this lineage and of the entire genus.

Environmental temperatures shape thermal physiology as well as diversification and genome-wide substitution rates in lizards.

Climatic conditions changing over time and space shape the evolution of organisms at multiple levels, including temperate lizards in the family Lacertidae. Here we reconstruct a dated phylogenetic tree of 262 lacertid species based on a supermatrix relying on novel phylogenomic datasets and fossil calibrations. Diversification of lacertids was accompanied by an increasing disparity among occupied bioclimatic niches, especially in the last 10 Ma, during a period of progressive global cooling. Temperate species also underwent a genome-wide slowdown in molecular substitution rates compared to tropical and desert-adapted lacertids. Evaporative water loss and preferred temperature are correlated with bioclimatic parameters, indicating physiological adaptations to climate. Tropical, but also some populations of cool-adapted species experience maximum temperatures close to their preferred temperatures. We hypothesize these species-specific physiological preferences may constitute a handicap to prevail under rapid global warming, and contribute to explaining local lizard extinctions in cool and humid climates.